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ABSTRACT 

We derive the anisotropic heat transport equation for rotating neutron stars 
and the thermal equihbrium condition in a relativistic rotating axisymmetric star 
is also derived through a simple variational argument. With a simple model of 
neutron star, we model the propagation of heat pulses resulting from transient 
energy releases inside the star. Such sudden energy release can occur in 
pulsars during glitches. Even in slow rotation limit (fi < 1 x lO^s"^), the 
results with rotational effects involved could be noticeably different from those 
obtained with a spherically symmetric metric in terms of the time scales and 
magnitudes of thermal afterglow. The effects of gravitational lensing and frame 
dragging on the X-ray light curve pulsations are also studied. The results we 
obtain indicate that the effect of rotation on the light curve modulation is small, 
and the spacetime curvature effect predominates. The metric components and 
rotational deformation on the stellar structure are calculated by using Hartle- 
Thorne formalism. We have applied our model to study the thermal 
response time scales of pulsars after glitches and find that the centrifu- 
gal force produced by rotation can substantially reduce the response 
time by a factor of three between non-rotating star and rotating star 
with ~ 900s^. The equation of state can also affect the duration of 
response. 

Subject headings: dense matter — stars: evolution — stars: interiors — stars: neutron- 
X-rays: stars 



^e-mail address: h9840329@hkusua.hku.hk 
^e-mail address: hrspksc@hkucc.hku.hk 



- 2 - 



1. Introduction 

The rotation-powered pulsars will spin down as they radiate. During the spinning-down 
epoch, some pulsars exhibit dramatic events call glitches, which is a sudden spinning-up of 
the star. Many observations have identified such phenomena (Krawczyk et al. 2003; Hobbs 
et al. 2002; Lyne, Shemar, & Smith 2000; Wang et al. 2000). There are numerous of mecha- 
nism have been proposed to explain the origin of pulsar glitches, while the superfluid-driven 
mechanism (Anderson & Itoh 1975; Alpar et al. 1984) and the starquake mechanism (Ru- 
dcrman 1969; Baym & Pines 1971) are most widely accepted. The starquake mechanism 
is based on the idea that a neutron star possesses a solid crust. As the star spins down, 
centrifugal force on the crust decreases and gravity pulls the crust to a less oblate equilib- 
rium shape. This change in stellar shape induces stress in the crust. However, the rigidity 
of the solid crust resists this stress and the shape remains more oblate than the equilibrium 
value. When the crust stresses reach a critical value, the crust cracks and the glitch energy 
is released in a small volume at the weak regions in the sohd crust. This locahzed energy 
releasing induces an uneven heating of the surface which corresponds to the 'spot' case. 
It should be noted that the starquake-driven glitches can occur anywhere within the crust 
where Coulomb lattice exists. 

The superfluid-driven mechanism (Anderson & Itoh 1975) suggests that the spin-up of 
the crust is produced by transferring angular momentum from a rotating superfluid to the 
more slowly rotating crust. Apart from spinning up, the mechanism also produces frictional 
heating and hence the local energy dissipation (Alpar et al. 1984; Shibazaki & Lamb 1988). 
Moreover, critical angular speed difference for spinning is inversely proportional to the dis- 
tance from the rotational axis, so only equatorial region shows concern to this process and 
hence a ring structure at the rotational equator is produced (Cheng et al. 1988). It should be 
noted that the superfluid-driven glitches can occur only in the inner crust where superfluid 
and Coulomb lattices coexist. 

Some other glitch mechanisms are proposed. For example. Link & Epstein (1996) have 
proposed a thermal glitch mechanism. In this model, a large increase in the vortex creep rate 
is induced by a temperature perturbation. As a result, the superfluid quickly loses angular 
momentum and delivers a spin-up torque to the crust. Carter, Langlois & Sedrakian (2000) 
have suggested that centrifugal buoyancy forces are the origin of pressure gradients sufficient 
to crack the crust and allow outward vortex motion. 

Serveral authors (Van Riper et al. 1991; Chong & Cheng 1994; Hirano et al. 1997; 
Cheng, Li, & Suen 1998; Tang & Cheng 2001) have suggested that the thermal evolution of 
a pulsar after glitches would provide a good method to determine the equation of state for 
neutron stars. The glitch energy was first assumed to be released in a spherical shell at a 
certain density inside the pulsar ( Van Riper et al. 1991; Chong & Cheng 1994; Hirano et 
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al. 1997), though this does not seem to be reahstic. They also show that the ghtches in the 
'shell' case cannot produce very significant observed results for young pulsars. 

Cheng, Li, & Suen (1998) argue that if a good fraction of glitch energy is released in a 
small volume, namely, the 'spot' case, then instead of heating up the entire stellar surface, 
even a small fraction of glitch energy can heat up a small area of the stellar surface drasti- 
cally. This would result in a periodic hard thermal X-ray pulse emission that should stand 
out clearly from the soft X-ray background. They suggest that by observing the time 
scales of thermal X-ray afterglows resulting from glitches, the equations of state 
of neutron steirs can be determined. However, the energy transport in the 'spot' case 
is clearly not spherically symmetric. They derive the general expression of the relativistic 
thermal transport and energy balance equations without assuming spherical symmetry and 
use these equations to study the evolution of the hot spot on the surface of neutron star 
after glitches. 

Tang & Cheng (2001) incorporate the relativistic light bending effect (Pechenick, Fta- 
clas, & Cohen 1983) and magnetic field effect (Page 1995) in their calculations. They show 
that these effects can significantly affect the intensity and the pulse shape of the transient 
X-rays resulting from glitches. 

Apart from the equation of state, glitch mechanism also affects the thermal responses. 
Larson & Link (2002) simulate the emergence of thermal wave at the stellar surface with two 
different glitch models and compare the results with the data from Chandra observations 
of thermal emission from the 2000 January glitch in the Vela pulsar (Helfand, Gotthelf & 
Halpern 2000). 

All the aforementioned studies have not considered the effects of rotating metric, which 
can affect the time scale, intensity and the pulse shape of the transient X-ray pulses. In the 
studies of 'spot' case (Cheng, Li, & Suen 1998; Tang & Cheng 2001), Schwarzschild metric is 
used and all the effects of rotating metric are ignored. For a complete analysis, development 
of an anisotropic transport equation is clearly needed and instead of Schwarzschild metric a 
rotating metric must be used. 

We organize the paper as follows. In section 2, we review the Hartle-Thorne formalism. 
In section 3, we study the effects of rotation on thermal equilibrium configuration. In section 
4, we describe the development of the anisotropic heat transport equation. In section 5, we 
describe a Markovian random walk method to simulate the thermal afterglows result from 
glitches. In section 6, we study the effects of the gravitational lensing and frame dragging on 
the X-ray light curves. In section 7, we discuss the physical interpretation of the numerical 
results. 
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2. Hcirtle-Thorne Formalism 

The general expression for the hne element of an axial symmetric spacetime is deter- 
mined by time-translational invariance and axial-rotational invariance: 

ds^ = e^^^'-'^Ut^ - e^^^'-'^^dr^ - r^e^^^^'^) [dO^ + sin^ ^ - u; (r, 6) dtf] (1) 

We have chosen units G = c = 1. 

To calculate the rotating stellar model as a perturbative expansion from a spherical star, 
we adopt the method proposed by Hartle (1967) and Hartle & Thorne (1968). We calculate 
the perturbations up to second order of the rotational frequency. The perturbed geometry 
of spacetime is described by: 



ds^ = e^"^''^ [1 + 2 {ho + h2P2 {cos 9))] dt 



[1 + 2 (2mo + maPa (cos &)) / {r - 2M)] ^ 



1 - 2M/r 

- [1 + 2 (i/2 - h2) P2 (cos 9)] [d9'^ + sin^ 9 {dcf) - udtf] 
+ 0(1]3) 



(2) 



where Vt is the rotational frequency of the star, P2 {cos9) = (3cos^^^ — l)/2 is the Legendre 
Polynomial of order 2 and hQ,h2,mQ, 7712, 1^2 are all functions of r that are proportional to 
Q'^ (Hartle 1967; Hartle & Thorne 1968). u {r) and M are the metric function and the 
gravitational mass of a non-rotating star respectively. With the chosen equation of state and 
the central density, they can be calculated by integrating the following set of equations 

— — = Anry (4) 
dr 

outward from the center with the boundary conditions: M (0) = 0; P (P) = 0; u (P) = 
i In ^1 — ^^^^^ j , where P is the radius of the non-rotating star. For a uniformly and rigidly 
rotating star, the contravariant components of four velocity are: 



1 + -r^ sin^ eo^e-'^'' -ho- /isPa (cos 9) 
2 



(6) 



(7) 
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u'^u^^O (8) 

where cu — — cu and cu defines the dragging frequency of local inertial frames. 

The star is centrifugally deformed by rotation. In the reference frame that is momen- 
tarily moving with the fluid, the energy-density distribution including the effect of rotation 
is: 

p + Ap^p+{p + P){p*+ p*P, (cos 6)) ^ (9) 

All the necessary rotational perturbation functions are calculated from the equations as 
described below. 

The angular velocity of the fluid relative to the local inertial frame is found by 
integrating the differential equation: 

1 d f ^ .doj\ , 4 (ij 



^-T r j-r + -- = (10) 
ar \ ar I r ar 



where 

3 (r) 



^ 2M 



(11) 



This equation can be integrated outward from the center of the star with the boundary 

dr 



conditions u — Uc and 4^ = 0. The constant Uc is chosen arbitrarily. When the surface is 



reached, one can determine the angular momentum J = (^)^_^ as well as the angular 
velocity fl = cu {R) + corresponding to the chosen cuc- 



Pq and mo can be calculated by integrating the following equations 

dmo ^dp 1 .2 ^ f duV 1 ^df _2 

— {p + P)po + —j r — --r—uj (12) 



dr dP ' " 12'' \dr ) 3 dr 



dpi _ mo(l + 87rr2p) 47r(p + P)r 



2 



dr (r - 2M)' (r - 2M) 



-Po 



Vl{r-2M)\dr) ?>dr\r-2M 

(13) 

outward from the center with boundary conditions Pq = mo = 0. Owing to the centrifugal 
force, a rotating star will have a larger mass than that of the non-rotating one for a given 
central density. The increment is given as AM = toq (R) + hg can be obtained from the 
algebraic relations 
Inside the star: 

ho^-pl + -r^e-^''0^ + hoc (14) 
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Outside the star: 



AM 



+ 



J2 



(15) 



r -2M r3 (r - 2M) 
hoc is a constant determined by demanding ho be continuous across the surface. 

The remaining four perturbation functions (z/2, ^2, ''^2,^2) t)e calculated from the 
following equations 



di'2 
dr 



2^^ h 
dr ^ \r dr 



3 dr G"' \dr ) 



(16) 



(iz/ 2r / (iz/ 
-2^ + 



dr r - 2M \ dr 

1 



27r (p + P) - 



M 



2z/2 



r{r-2M) \dr 



du 



+ 



di/ 



-r — 



dr 2(r-2M) V^?^ 



dz/ 



-1 



dc<j\ ^ 



dr / 



du 

—r + 



dr 2(r-2M) \dr 



du 



2-2 ■ 

r 00 J — 
dr 



m2 = (r - 2M) 



—ho r ' 1 . I -''-^ 



dr y 6 V dr 



* U 2 

P2 = -^2 - 3^ e 



with boundary conditions U2 — and /i2 = at the origin. 



(17) 
(18) 
(19) 



3. Thermal Equilibrium Configuration 

Consider a particle in the star with energy E, the energy of this particle as measured 
by a distant observer E^o can be related to its locally measured value E by E^o = Next, 
we consider the energy transport from one region to another in an isolated star without any 
change in the density of particles. For a local observer, the fundamental thermodynamics 
relation dE — TdS is always vahd, where dS is the change in entropy per baryon and T is 
the locally measured temperature. 

By following the scheme of Zel dvoich and Novikov (1971), the equilibrium condition 
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can be obtained by extremizing the entropy of the system with the constraint that the total 
energy is conserved: 



{/■ 



SndV + A / E^ndV = 



SS + A^]ndV = 
) 

1 + 1 ndV = 



(20) 



where n is baryon number density, dV is proper volume element and A is a Lagrange multi- 
plier. Consequently, the thermal equilibrium condition of a relativistic axisymmetric star is 
formulated as: 

T 

— = A~^ = constant (21) 

This result, which obtained from a relatively simple variation argument, is consistent with 
the one derived more vigorously by Miralles, Van Riper, & Lattimer (1993). Also, there is an 
underlying assumption that the neutron star is in rigid rotation. If it were not, the internal 
friction will produce heat and our variational argument is no longer valid. 

The aforementioned isothermal approximation is only appropriate in the interior. There 
is a thin atmosphere surrounds the isothermal core which sustains an appreciable tempera- 
ture gradient. For spherically symmetric cases, model atmospheres have already been calcu- 
lated by Gudmundsson, Pethick, & Epstein (1983), which give the surface temperature Tg 
as a function of the temperature at the base of the atmosphere and the surface gravity gg. 
The base of the atmosphere is chosen at a particular density ph- For our study, it is chosen 
at pb = 10^° g cm~^: 

(ji4 s 0.455 
^ ) (22) 

We shall assume this relation is valid in the slow rotation case with a suitable choice of 
an effective local surface gravity gg iff) that take the centrifugal deformation into account. 
Consequently, the surface temperature Tg {Ti,, 0) depends on the polar angle. 



4. General Relativistic Anisotropic Heat Transport Equation 



With an intention to model the most general case of energy transport inside a rotat- 
ing neutron star, we derive the heat transport equation without making the assumption of 



-8- 



spherical symmetry in energy transport, stellar structure as well as the spacetime. 

To derive the equation, we assume that there are no other entropy generating mecha- 
nisms besides diffusion and the ffuid motion inside the star due to the thermal effects is neg- 
ligible. The energy-momentum tensor inside a star consist of perfect fluid which allow heat 
flow can be written as T'^'^ = {p + P) u^u'^ — P g^'^ q'^ ^ where g^^ denotes the inverse 
components of equation (1), is the 4- velocity of the fluid flow which has been described 
in section 2, p and P are the total energy density and pressure measured in the rest frame 
of the fluid(Tolman 1934). The heat flow is given 

where = Ua-pu^ is the 4-acceleration and K is the thermal conductivity(Tolman 1934). 

No matter how complicated the transport equation is, the basic principle is that the 
underlying physics is nothing but the law of conservation of energy. In relativity, the energy- 
momentum tensor embodies a compact description of energy and momentum. To be more 
speciflc, the transport equation can be derived by the conservation of energy-momentum 
tensor, namely. 



where 



t: 



ta 



Tf + T. 



te 



ntcf) 







(23) 



T: 



+ 
+ 



t: 



dP 
dp 



dp 

at 



u u 



dtdcj) 



tr 



8f 



t: 



rptrj) 



Of 
09 
df4> 



(24) 
(25) 
(26) 
(27) 



where /* = T is the energy flow per unit area parallel to i direction. F^^ are the Christoffel 



symbols. || in is the rate of change of energy density measured by distant observer, 
which can be expressed in the locally measured quantity, namely, It depends on the 
processes under consideration. If only heat conduction and neutrino emission are considered, 
we have |^ = -\-Qi,, where Cy and Qi, are the heat capacity and the neutrino emissivity 
measured in the local frame respectively. Expressing it in the coordinate time t, we have 
|2 — Cy^ + ^. Using this relation and equation (23), we obtain the heat diffusion equation. 







dP 
dp 



at w 
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+ 2u^K 
df 



+ 
+ 
+ 



(,«_,.V)^ + (,'*-«V)|^ 



Or 



df df 



oe 

dX 
dr 
dX 
89 

+ 2r* A'' + 2rUV 



f 



(28) 



In the non-rotational limit, the spacetime is described by the diagonal Schwarzschild 
metric, the above equation recovers the one obtained by Cheng, Li, & Suen (1998). It should 
also be noted that this equation recovers the well-known Newtonian case. 



5. Simulations of Thermal Afterglows result from Glitches 

The anisotropic relativistic heat transport equation, i.e. equation (28), can 
be rearranged into a single variable partial differential equation (cf. Appendix). 

We employ a quantized Monte Carlo technique for the simulation of heat transport inside a 
rotating neutron star. Dynamic stochastic processes are simulated by using rate coefficients 
(i.e. diffusion coefficient D, drift coefficient v, annihilation coefficient /i and local depletion 
S). These coefficients are determined via rearranging the rate equation (28). In order 
to perform the simulations, some assumptions must be made. First, we assume that the 
quantities D^, Vi and /li depend only on position where i denotes the motion of energy carriers 
in r, 6 and 0. Much more general dependences result in non-linear behavior. Second, if the 
grid sizes are fine enough, the energy carriers transit from the surrounding grid points in 
the previous time step with equal probability. Third, we neglect the thermal effects on the 
stellar structure. 

There are two terms of mixed derivative and second time derivative present in the 
equation. There is no easy way to incorporate them in a random walk approach. However, 
their coefficients only contain the factors of (g'** — u*u*) and {g*'^ — u^u'^) . Since the zero order 
terms are cancelled in the subtractions, only second order perturbation terms are remained. 
Up to the highest rotational frequency we have considered (1 x lO^s""^), these coefficients 
are still negligible in comparison with the smallest term we have taken into account. So it 
is safe to drop it in the simulations. 

With these assumptions, the heat transport inside a rotating neutron star can now 
be simulated by a straightforward Markovian random walk (Mackeown 1997): (1) Specify 
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the initial condition T (r, 0) and quantize it by a large number, M, of "walkers" whose 
starting positions are selected in accordance with T (f, 0) and their random walk is modeled. 
(2) In each time step At, a weight is attached to each "walker" according to the survival 
probability fiAt by means of survival biasing (Mackeown 1997). Then, it is displaced by 
Axi = ±V2-DjAt + ViAt in each degree of freedom, where the sign of the step is chosen 
with equal probability on the basis of a uniform random variate. (3) Incorporate the local 
depletion SAt. (4) Repeat steps 2 and 3 M times for each "walker". (5) Repeat step 4 
for N times. (6) Plot the distribution of the "walkers" and this should approximate the 
solution T (r, t) at t — NAt. Note the finite value of "walkers" M introduces an unavoidable 
statistical fluctuation. Also, the flnite time step used will introduce a truncation error. 

Since we intent to investigate the effect of rotation on the heat transport, we adopt a 
simple model of neutron star. We employ the method described in section 2 to calculate the 
rotational stellar structure. The stellar structure is determined by the equation of state. We 
use the equation of state of neutron matter from Pandharipande (1971) in our calculation. 
We have considered the contribution of protons, neutrons and electrons to the heat capacity 
(Maxwell 1979). But the most important term is the contribution of electrons. We only 
consider the neutrino emissivities due to neutron-neutron bremsstrahlung, proton-neutron 
bremsstrahlung and the modifled URCA process (Maxwell 1979). We adopt the analytic 
formulae of thermal conductivity presented by Flower & Itoh (1981). 

We calculate temperature distributions for 'ring' case and 'spot' case. For both cases, 
heat inputs are deposited at the depth of the crust where p ~ 10^'^ g cm^'^. The ratio 
R/M ~ 8 corresponds to M ~ IMq and i? ~ 10 km. The core temperature is taken 
to be lO^K. With AE = lO'^^ergs deposited around the equator (i.e. 9 = 90°), figure 1 
illustrates the polar angular surface temperature distribution of a rotational 'ring' case with 
Q — 7.6x 10^s~^. For the 'spot' case, we choose AE — lO^^ergs to be released at 9 — (f) — 90°. 
Figure 2 illustrates the azimuthal profile of surface temperature at the rotational equator of 
the hot spot for a rotational case with Q = 7.6 x 10^s~^. Figure 3 illustrates the evolution 
of thermal X-ray fiux of 'spot' case. Three cases with different rotational frequency are 
compared. With the mass increased by ~ 33% and radius decreased by ~ 1%, the R/M 
ratio is reduced to ~ 6 which indicates the gravitational effect is enhanced. We re-calculate 
the 'spot' case with this R/M ratio and the results for three cases with different rotational 
frequency are shown in figure 4. The time scale and the magnitude of the afterglows for the 
non-rotational cases are consistent with those obtained by pervious authors (Tang & Cheng 
2001). 

The heat propagation time r is approximately inversely proportional to the diffusion 
coefficient Dr, which is proportional to thermal conductivity K and inversely proportional 
to heat capacity Cy. In the region we are interested, the heat capacity is proportional to 
temperature, while the thermal conductivity is inversely proportional to the temperature. 
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Hence, r is roughly proportional to the temperature squared. It can be seen that heat 
propagation becomes much slower as the temperature becomes higher. For this reason, the 
duration of afterglow for the 'ring' case is shorter than that for the 'spot' case for the same 
heat input. Comparing the non-rotational flux curve of 'spot' case with R/M ~ 8 in figure 3 
and that with R/A4 ~ 6 in figure 4, wc find that the afterglow take a longer time in the model 
with R/M ~ 6 because a larger mass underneath the location of energy release that hinder 
the propagation of heat pulse to the surface. Comparing the rotational cases with the non- 
rotational one in figure 3 & 4, we find that the duration of thermal afterglow is shortened 
when the effect of rotation is introduced. To further investigate this effect, we carry out 
more simulations of 'spot' case with different rotational frequency. Figure 5 illustrates the 
fractional decrease in the duration of thermal afterglow as a function of rotational frequency 
for R/M ~ 8 and R/M ~ 6. 

The rotational effects can be originated from the centrifugal forces in rotating neutron 
stars. When the energy Ai? is released in a localized region at the equator, it can be 
viewed as an equivalent mass ~ IS.E /(? that corotates with the star. This leads to a natural 
interpretation that centrifugal forces prompt a faster and larger thermal response. In figure 
5, it should be noticed that the fractional decrease in the duration is larger in the case of 
i?/M ~ 6 than that in the case of it!/M ~ 8 for a given fi. This indicates the other effect of 
rotation on the heat transfer. Apart from prompting a faster heat propagation, centrifugal 
forces also deform the star. As pointed out by some ellipticity studies (Abramowicz 1990; 
Chandrasekhar & Miller 1974), the denser a neutron star, the less eccentric it is. The star 
with R/M ~ 8 is more eccentric than the one with R/M ~ 6. Since the heat inputs are 
deposited at the depth where p ~ lO^^g cm~^ at the equator in both models, the heat pulse 
needs to propagate through a larger distance in a more eccentric star in order to reach the 
surface and give the afterglow. 

Nevertheless, the behavior of centrifugal forces in general relativity is fundamentally 
different from that in Newtonian physics (Freire & Costa 1999; Abramowicz, Carter, & 
Lasota 1988; Abramowicz & Prasanna 1990; Abramowicz, Miller, & Stuchlik 1993). One 
of the questions behind our interpretation is whether there can be inversion of centrifugal 
forces exist inside neutron stars. Freire & Costa (1999) show that most realistic equations 
of state do not allow the existence of such inversion. Without the bother of centrifugal force 
inversion, we suggest that rotation introduces two effects on the thermal afterglow. First, 
it prompts a faster and larger thermal response. Second, it increases the distance that heat 
need to travel through in order to give an afterglow on the surface. The resultant effect on 
the duration and the magnitude of the thermal afterglow depends on the interplay between 
these two factors. 



6. Thermal X-ray Light Curves 



For the 'spot' case, the locally released energy would modulate the X-ray pulse shape 
by heating a portion of crust so that more thermal X-rays are emitted at a particular phase. 
The modulation of X-ray pulse will last until the surface temperature equilibrated. We have 
investigated the effects of spacetime curvature and rotation on the thermal X-ray profile. 



6.1. Gravitational Lensing Effect 

Since the surface gravity of a neutron star is tremendous, the effect of gravity on the 
trajectory of emitted photons must be taken into account. In this section, we choose our 
coordinates so that the observer is on the positive 2;-axis at r = Tq where Tq —>■ oo (see figure 
6). The following scheme of calculating the X-ray light curves resulting from light bending 
is adapted from Pechenick et al. (1983). The surface of the star is described by angular 
spherical coordinates, 9 and 0, where 9 is measured from the z-axis we have just defined. 
When the photon is being emitted from the surface at an angle S as illustrated in figure 6, 
it will be deflected by the gravitational field. It will seem to the observer that it is emitted 
at angle 9' from the z-axis. Hence, ^ is a function of 9': 



T 



- (1 - 2u) u 



du (29) 



where h = rQ9' and u = M/r. b is the impact parameter of the photon. 

If the star has the ratio R/M > 3, then a photon emitted from the surface and reach 

the observer must have an impact parameter b < b„iax, where bmax = R{l ~ ^) ^ The 
maximum possible value of 9 occurs when b = bmax- We consider a hot spot of angular 
radius a centered at 6* = 9o. A function h{9; a, 9q) is then defined as the range of included 
in the "one-dimensional shce" at 9 of the hot spot. II 9o + a < 9 max < 180° and 9q — a>Q, 
then h{9] a, ^o) is defined as: 

f 2cos-^ /" cosa-coseocosA - a < 9 < 9o + a; 

h{9]a,9Q)^< \ s'-i^osine ; _ _ u , 

[ for ^ outside the range 9Q±a 

When the requirement that the photons under consideration must reach the observer is 
imposed (not all the photons emitted at (^,0) can reach the observer), 5 is also a function 
of 9' . Let X — and Xmax — we have 
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We will come back to this requirement when we discuss the rotational effect on the pulse 
shape. 

In order to obtain light curves, has to be obtained as a function of time. If j3 is the 
angle between the axis of rotation of the star and the line joining the center of the hot spot 
and the center of the star, and 7 is the angle between the axis of rotation and the 2;-axis (see 
figure 7), then 

cos 00 = sin (P) sin (7) cos {Vlt) + cos {/3) cos (7) (32) 

where Q is the rotational frequency of the star. We are now able to calculate the relative 
brightness: 

A {9o; /, M/R, a) = (l - — J ( — ) f {S (x)) h {x; a, Oo) x dx (33) 

where f {5) — 1 for isotropic emission, / [5) — cos 5 for enhanced emission and / {5) — sin 5 
for suppressed emission. The relationship between A and is plotted in figure 8. 



6.2. Rotational Effect 

Apart from getting bent due to spacetime curvature, the trajectory of a photon emitted 
in a general direction from a point rg will be dragged away from its original direction of 
emission when the rotation of the relativistic star is taken into consideration (Kapoor & 
Datta 1985). In this section, the axis of rotation is taken to be z-axis. The net angle of 
defiection in the trajectory will be given by: 



e'^-^ [(1 + ujqe) - qle^^ - e^^] ' 

with Te and Tq denote the point of emission and the observer's location respectively, e*^, 
and are the exterior rotating metric components: 

e''^ = sin^ d (36) 
The dragging frequency uj has an analytic exterior solution: 

a;(r) = ^ (37) 
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ge in the integrand of equation (34) is defined as the impact parameter of the photon with 
rotation taken into consideration: 



1 + e^-*^ {e-^-^u {n-u;) + Q sin 6) 



(38) 



In our case, the source (hot spot) is located at the rotational equator, 6 represents the 
azimuthal angle at which the photon is emitted with respect to the normal vector of the 
surface as seen in the local rest frame of the star. Wc choose the convention that S = for a 
radial outgoing photon, 5 < for a tangentially forward photon, and (5 > for a tangentially 
backward photon. 

As mentioned in the section 6.2, only the photons emitted at a particular value of 5 
would be received by the distant observer. The value of qe has been obtained by imposing 
this requirement. Because of the effect of frame dragging around the neutron star, the 
situation will be different from that in previous section. 

It should be noted that the impact parameter is not symmetric in ±S. In principle, 
this asymmetry will manifest itself in the final pulse shape: 

5neu,^S + (l>o{S) (39) 

Figure 9 schematically illustrates the deflection of a photon trajectory. As a consequence of 
rotation, |0o(^)| |0o(-'5)|- 

This asymmetry does not occur in spherically symmetric spacetime, where although the 
spacetime curvature will still deflect the trajectory of photon. This can be easily seen if one 
approximates the corresponding Schwarzschild expression of equation (34) as follows: 



— 'r 



sin^ S {2R - 3R^) 
^'^ 12 {R- 2RS) 



(40) 



where the superscript S refers to the Schwarzschild case. Here 6 = 90° for the reason of 
symmetry. R^ denotes Schwarzschild radius and q^ — {e'^~'^ sin 8) It can be seen that 

We define the corresponding 5^^^^ of equation (39) as 

€ew-5 + ^U5) (41) 

comparison between 5^^^ and 5^,^^ is made in figure 10. Even in the fastest rotational case we 
have considered (Q = 1 x lO^s"-*^) with a relatively small gravitational effect {R/M ~ 8), the 
comparison indicates that the effect of rotation is small and the effect of spacetime curvature 
predominates. 
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6.3. X-ray Light Curves 



Wc can generate the light curves in soft X-ray regime for different 'spot' cases. According 
to Pechenick et. al. (1983), the total energy flux observed is: 



where Iq is the energy flux from each cell at the surface including the factor of enhanced 
emission. We adopt 5new instead of 5^^^ to incorporate the asymmetry due to rotation. ^ 
denoted the summation of the contribution from each cell. The comparisons of light curves 
with different rotational frequency are made in figure 11 and figure 12 with R/M ~ 8 and 
i?/M ~ 6 respectively. We have not found any asymmetry in the pulse profile for the rota- 
tional cases. 



In this paper, we have developed an anisotropic heat transport equation. This equation 
is then employed to simulate thermal afterglows resulting from pulsar glitches. We examine 
the effects of rotation on the thermal responses. Although rotation does not have significant 
effect on the rate of standard cooling (Miralles, Van Riper, & Lattimer 1993), noticeable 
changes on duration and intensity of thermal afterglows are found in 'spot' cases. In com- 
parison with static cases, afterglows with shorter duration and larger intensity are found in 
rotating stars. These are not unexpected when the effects of centrifugal forces on the stellar 
structure as well as heat transport are realized. We suggest that rotation prompts a faster 
and larger response. In equation Al(cf. Appendix), the radial diffusion coefficient 
is given by Dr = which contains the factors of rotating metric and hence is 
a function of rotational frequency. Dj. increases with rotational frequency (see 
table 1) and prompts a faster response. On the other hand, the rotation also increases 
the distance that heat pulses need to travel through in order to give an afterglow on the 
surface. The resultant effect on the duration and the magnitude of the thermal afterglow 
depends on the interplay between these two factors. We have also generated the thermal 
X-ray fight curves for the 'spot' cases. The effects of gravitational lensing and frame dragging 
are fully taken into account. Apart from the peak time and the peak intensity, we have not 
found any significant effect of rotation on the morphology of the pulse profile. 

Hirano et al. (1997) characterize the thermal afterglow by the fractional increase of 




(42) 



7. 



Conclusion and Discussion 
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surface temperature at the peak and the peak time. They also propose a framework to set 
constraint on the equation of state by using these two parameters. Since the peak time is 
larger and the amplitude is smaller for a stiff star than a soft star. Without taking the effects 
of rotation into consideration, we may underestimate the stiffness of the equation of state. 

One of the practical questions is the detectability of these soft X-ray transients which 
depends on the sensitivity of the detector. For the 'spot' cases considered in this paper, the 
glitch events are able to be detected by the state of the art X-ray sateUites such as XMM- 
Newton. Once the observational data of these thermal afterglows are obtained, these can 
be used in putting constraints on the equation of state of neutron stars as well as the glitch 
models by the method of periodic analysis (Andersen & Ogelman 1997). 

In our calculations, we have neglected interstellar absorption, magnetospheric effects, 
magnetic field effect (Page 1995) and the possibility of uplifting and local expansion of mat- 
ter caused by the heat deposition (Eichler & Cheng 1989). We would like to remark that 
we have ignored the contribution of impurity scattering in the thermal conduc- 
tivity, which may be important when T < 10^ (Yakovlev &z Urpin 1980). We have 
also assumed that there is no other glitch events occur during the afterglow is evolving. If 
the glitch events occur too frequently, it may result in a pile-up of pulses as well as a long 
term variation of the total thermal radiation which eventually reduces the detectability of 
the thermal afterglow (Li 1997). Moreover, we keep a constant rotational frequency during 
the afterglow is evolving. However, glitch recovers exponentially and hence Q should varies 
with time. Also, the effects should be more important for fast rotation, a fully numerical 
scheme is needed to calculate the metric in fast rotation case. For further study, all these 
effects have to be taken into account. 



The anisotropic relativistic heat transport equation (i.e. equation (28)) can be rear- 
ranged into a standard form of partial differential equation as follows: 



A. Appendix 




Q2rp Q2rp Q2rp g2rp g2rp Qrp QJ. QJ, Qj, 



Ci = 2u*K (y" - u^u* 



) 



(A2) 



C2 = u'Kg 



rr 



(A3) 
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C3 = u'Kg'' 
C5 = K [3m* {g'^ - mV) + (^** - 



ax 



^(u'u'-g'')+u'u' 
dp 



(A4) 
(A5) 
(A6) 

(A7) 



dr 



dr 



dX 



Cr = K—g-+u'—g-+u'K4-+ ( 3r*, + ^ + r% + tJ, ) u'Kg-+2rl,u'^Kg--u'Kg-a 



dr 



de 



de 



de 



(A8) 



^8 = i^^/'^ + u'-^-^g'' + M*i^ ^ + ( Sr*, + ^ + C + r^, M^i^/*^ - u^Rg'^'ae (A9) 



10 



Oil — — — 



(AlO) 
(All) 
(A12) 
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Fig. 1. — Surface temperature of a 'ring' case as a function of polar angle for a neutron star 
with RjM ~ 8, = 7.6 X lO^s-i, = WK, AE = W^^eig at Pgiuch = lO^^g cm-^, 9 = 90°. 
The upper panel denotes days elapsed after the heat deposition. 
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Fig. 2. — Surface temperature of a 'spot' case at the equator as a function of azimuthal 
angle for a neutron star with R/M ~ 8, Q = 7.6 x = lO^K, /^E = lO^^g^g at 

Pgiitch — lO^^g cm~^, — (f) — 90°. The upper panel denotes days elapsed after the heat 
deposition 
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Fig. 3. — The evolution curves of thermal X-ray flux for a neutron star with R/M ~ 8, 
= lO^K, AE = lO^^e^g pg^■^^f^ = lO^^g ^^^-3^ q ^ ^ ^ qqo ^ Three cases 

with difi^erent rotational frequency are compared. 




Fig. 4. — Same as figure 3 but with R/M 6. 
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Fig. 5. — Fractional decrease in the duration of thermal afterglow of 'spot' case as a function 
of rotational frequency for R/M ~ 8 and R/M ~ 6. 




Fig. 6.— 



The geometry used to determine the brightness of a hot spot. 
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Fig. 7. — The relation of the orientation of the hot spot and the axis of rotation to the 
observer's hne of sight. 
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Fig. 9. — Schematic illustration of the effect of frame dragging on the photon trajectory. 
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Fig. 10. — Effect of frame dragging on the photon trajectory. illustrates the deflec- 
tion in Schwarzschild spacetime and S'^, S" illustrate the deflection in tangentially backward 
direction and tangentially forward direction at the rotational equator of a neutron star re- 
spectively. 
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Fig. 11. — Thermal X-ray pulse profiles at the peak time for a neutron star with R/M ~ 8, 
Te = lO'^K, = lO^^erg at pgutch = lO^^g cm-^ 6* = = 90° in an orientation with a = 3°, 
f (S) = cos 5, /3 = 7 = 90° for 'spot' case. Three cases with different rotational frequency 
are compared. The numbers associated with the arrows indicate the number of day after the 
glitch 
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Fig. 12. — Same as figure 11 but with R/M ~ 6. 
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i?/M = 6 




AD J. cm^s ^ 




ADj. crri^s ^ 


152 


14.1 


200 


28.6 


305 


56.3 


300 


60.3 


500 


152.2 


513 


183.0 


579 


201.3 


586 


238.0 


762 


344.6 


760 


395.8 


839 


416.1 


802 


418.2 



Table 1: The deviation of radial diffusion coefficient from tfie non-rotational value (i.e. 
ADj. — Dr^Q^o — -Dr,n=o) at the location of energy deposition for 'spot' cases (i.e. at 
Pgiitch = lO^^g cm-^ e = = 90° with AE = lO^^erg, = lO^K) at i = s as a func- 
tion of rotational frequency. 



